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1 Introduction 



One-dimensional models with inverse square interactions have attracted considerable interest 
in recent years. For the Calogero Sutherland (CS) models |T]-^| describing particles moving on 
a continuous line the many-body ground state wave function is of Jastrow type and excitations 
can be written as a product of this pair product wave function and certain polyniomials in 
the coordinates. Based on this observation the eigenvalues of the Hamiltonian can be found 
by means of an asymptotic Bethe Ansatz (ABA) solution 0. Finite-size scaling analysis of 
the excitation spectrum and predictions of conformal field theory (CFT) have been used to 
study the critical behaviour of the CS model, leading to the identification of the universality 
class of the model with periodic boundary conditions as Luttinger liquid, i.e. a Gaussian model 
with central charge c = 1 j§,|5[]. An interesting property of th CS models is that the compact 
form of the eigenstates allows for an explicit calculation of certain correlation functions [0 0. 
thus allowing to compare the asymptotic predictions of CFT to exact expressions derived in a 
microscopic model. 

In addition to the models with periodic boundary conditions there exists a class of CS 
models lacking translational invariance, in particular the model of -BCjy-type which is invariant 
under the action of the Weyl group of type M. Again the ground state can be written 



in a compact form of Jastrow type, and the spectrum of excitations can be found [10,111 



An analysis of the finite size corrections to the energies shows that the spectrum acquires 
contributions due to the 'boundaries' of the system and the low-energy critical behaviour is 
described by a c = 1 boundary CFT [0. Again, the existence of a 'simple' expression for the 
ground state wave function opens the possibility to compare the predictions of boundary CFT 
for the asymptotic behaviour of correlation functions with exact results. These predictions are 
of great interest at present due to the possibility to extract observable properties of quantum 
impurity systems from finite size spectra (see e.g. fll3|-|T5|| ) . 

In the present paper we shall consider the -BCV-type CS model and compute matrix el- 
ements which allow to compare the exponents associated with the Anderson's 'orthogonality 
catastrophy', i.e. the dependence of the overlap between ground states for different boundary 
conditions on the system size, and the asymptotic behaviour of (Friedel) density oscillations 
due to the existence of the boundary with the corresponding expressions for the Luttinger 
liquid ||H2g . 

First we will give a brief review of the properties of the CS model of -BCjv-type as well as the 
predictions of boundary CFT relevant to the results of this paper. A form of the Hamiltonian 



in a finite geometry especially convenient for our studies is [VI 
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Here A, /x, i/ are positive coupling constants. The particles move on the interval < qj < L. The 
two particle interaction terms consist of the usual inverse square interaction of particles moving 
on one half of the circle with circumfence 2L plus a term from the interaction of the particle 
at qj with the mirror image — q k = 2L — q k of the particle at g&. The last two terms can be 
regarded as impurity potentials situated at the edges q = (q = L) of the system with strength 
given by /i {y). Eigenvalues and eigenstates of (|l|) have been determined explicitely [10, 11 
the many-particle ground state wavefunction is M 
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As for the periodic CS models the spectrum can be reproduced exactly by means of the asymp- 
totic Bethe Ansatz method [Kj. Expanding the ground state energy in inverse powers of the 



system size the following finite size scaling form is found [12 



E8> = L e <°> + 2/ + ^W) 2 -^fA (3) 
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where e*- - 1 and / are the bulk energy density and the boundary energy in the thermodynamic 
limit for fixed particle density n = N/L, respectively. Vp = 2irXn is the Fermi velocity of the 
elementary excitations. Similarly, the energy of an excited state with AN additional particles 
and N p h > particle-hole excitations near the Fermi point is, to leading order in l/L, 



E - < - „<°>AiV = ^ Q (AN + AN„f - ± (AN„f\ + 



^N ph . (4) 



Here we have absorbed a term AJV into the definition of the energy the system with /j,^ = 
Lde^ /dN being the chemical potential. 

The expressions (§), (f|) should be compared with the corresponding prediction of CFT for 
models with free boundary conditions, namely 

B W = ie (») + 2/-^c + ^x (5) 



with the Virasoro central charge c appearing in the universal amplitude of the 1/L-term and 
the critical exponent x of the operator generating this state. It is well known that the long- 
range nature of the interactions in CS models gives rise to non-universal 1/L-contributions to 
the ground state energy of these systems |j5, 12|. The A-dependence in the last term in (|||) 



is believed to be a consequence of this effect, thus yealding an incorrect value for the central 
charge. 

The other contribution of order 1/L to the ground state energy is a consequence of the scat- 
tering due to the free boundary and the impurity potentials. If one applies boundary CFT to 
obtain the surface critical exponents controlling the asymptotic behaviour of correlation func- 
tions one has to distinguish operators connecting states corresponding to different boundary 
conditions and operators inducing a change of particle number or particle-hole excitations in 
the ground state corresponding to a given boundary condition |]i4|,|l5|j. For the latter case the 
phase shift AN b should be absorbed into the the change of the number of particles 

Ef = 4 0) - ^ m b ) 2 , AN = AN + AN b (6) 

to restore particle hole symmetry of the finite size spectrum (|j). The resulting scaling dimension 
of an operator <fi corresponding to this situation is 

with integer AN. 

To obtain the conformal dimension of boundary condition changing operators finite size 
energies corresponding to states subject to different boundary conditions have to be compared 



Tj.[Ij||17] . Consequently, only one of the two phase shifts AN b , ANy arising in these expressions 
can be absorbed into a shift of the particle number as in (|6]). For the operator ifj bb > connecting 
the gound states corresponding to different boundary conditions this leads to an operator 



dimension 



x(ipw) = - (AN b - AN b ,) 2 . (8) 



This exponent determines the orthogonality exponent (0 b \0 b ') oc L~ x and the related X-ray 
edge singularity arising from a sudden change of the boundary potential |l],[l£||17|,|21j. 

The asymptotic behaviour of correlation functions in the bulk is still determined by the con- 
formal dimension of the corresponding operator. In fact, an n-point function of the semiinfinite 
system is subject to different boundary conditions but obeys the same differential equation as 
the 2n-point function including the mirror positions in the system without boundary ||22|| . 



Hence, the critical exponent for the asymptotic behaviour of the single particle density (p(q)) 



due to the presence of the boundary can be obtained from the density density correlation func- 
tion for the corresponding system without a boundary, namely the periodic CS model. The 
most dominant term is due to backscattering processes with momentum ±2&f = 27m and 
decays asymptotically as f|: 

cos(2k F q) 

(p{q))-n -qTJx—- ( 9 ) 

2 Overlap integrals 

To compute the overlap integral between ground states of ([!]) corresponding to different values 
of the boundary field stengths /i, u we will make use their representation in terms of so called 
Selberg correlation integrals (see Hj and references therein) 
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For m = the integrals can be evaluated with result expressed in terms of Gamma-functions 

o n x ^ _ rr r d + xmrjXi + Xjj - i)/2)r(A 2 + A(j - i)/2) 

Snfi{Xu A2 ' A) " lj r(l + A/2)r(A 1 + A 2 + A(n + j-2)/2) ' {U) 
We first note that the normalization integral of (|2|) 

Mx,w= \Jl£ d 9ij \* { o\qi,...,qNlH,")\ 2 (13) 
is of the form (]T0|) : Substituting ti = \ (l — cos^qij we obtain 



i2A 



l<j<fc<n 
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Now the overlap between two states (0) for different sets of boundary fields (/x, z/) and (//, z/) 
can be expressed as the ratio of Selberg correlation integrals 

i, ,, « .2 GW^ + ^ + x )/ 2 ' (" + 1 + ^ 2A )) 2 n ^ 
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Using ( |T2Tj we obtain 

1 JV ~ 1 1 1 

In | (//, v\n\ z/) | = - (/O, /A 2 + X fi + /("> ^' 2 + 

3=0 

-/ + */,// + z/, 1 + A(iV + j - 1))) (16) 

with 

/(,, »> ,) = h r r ;^f x) + r y4 « o. - (-i + ^ o. + ,/ - d + ■ ■ ■) . (it) 

For large system size L and fixed density p = iV/L the overlap integral (|16|) is determined by 
the logarithmic divergence of the first two contributions to the sum at j = N — 1 giving 

In u\n', u')\ oc -L ((/i - //) 2 + (i/ - z/) 2 ) InL (18) 

which is in perfect agreement with (|3|) and the prediction (§) of boundary CFT (note that the 
operator connecting v) and |//, u') in this situation is a product of two boundary changing 
operators ip^' an d acting at the boundary at q = and q = L, respectively [FT]). 



3 Friedel oscillations 

We consider now a spatial dependence for single particle density 



For integer values A the function (p(q)) can be expressed in terms of the Selberg integrals fllQ|) 

L S Nfi (// + 2' 2A J 

where m = 2 A and X\ — . . . — x m — x — sin 2 (7rg/2L). Eq. (|20"D can be rewritten as 

£ <SjV,o(/4 + 2'^+ 2' 2A ) 

1 11 

x 2 F[ X \-N + 1, -(/i + i/ + m) + iV - 2, -(// - - + m), x x = . . . = x m ), (21) 

where is a generalized hypergeometric function of m variables (see 0). For finite systems 
this expression for the single particle density can be evaluated by using the fact that 2 -^i f° r 
equal arguments can be written in terns of Jack symmetric polynomials p3 |. In Fig. [I] we have 
plotted ( PH) in this representation for a system of N = 10 particles. 



In the thermodynamic limit we derive the asymptotic behavior of (p{q)) for 1 <C q <C iV 
using the integral representation |7| 

2 F 1 (2/A) (a,Ai + A(m-l)/2,A 1 + A 2 + A(m-l),xi = . . . = x m ) 

I [ / dt t (l - xUY a Dx^xih, ...,t m ) . (22) 



Omitting the ^-independent factor we find from ( pT| ) 

w^oc^-.rn U *{*£ffi} ^J nto-^r. ^ 



where n = iV — 1. Near the boundary g< 1 we obtain 

(/>(?)) oc . (24) 

For 1 C g < Af we have x = sin 2 (7rg/2L) — > g 2 /n 2 , g = kpq/2. After rescaling t; — > nt; 
and using the identity lim^oo (1 + y/n) n = expy we obtain from 
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In order to calculate this divergent integral we analytically continue it to imaginary g and 
obtain after rescaling t —*t/g 
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j [ / c^exp-^+^tf- 4 ^ )l[\t j -t k f\ (26) 
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Now Eq. (|26|) can be easy evaluated in the limit q ^> L/N which corresponds to g ^> 1. 
Considering the case m = 2 first we find after an integration near the extremum point tj = 1 

cos2k F q 

In the general case the main oscillating term is obtained by integration ( p6| ) over an region where 
only two variables are near extremum points. Then the integration over these two variables 
gives rise to a factor 

cos2kpq 

gl/A+l 

while the other integrals of type Yl j e ~ gt f{t) contribute a factor of order (^) m_2 . Substituting 
to Eq. fl26|) we reproduce the result ( p7| ) for any value of the coupling constant A. Again, this 
result is in complete agreement with the boundary CFT prediction (0). 



4 Free fermionic case 

Note that in the free fermionic case A = 1 the Friedel oscillations can be investigated in more 
detail due to a possibility to express the Selberg integral S n ^ in terms of Appell's hypergeometric 
function F 4 

S n ,2{M, A2, A; xi, X2) 

= (-l^Sn^ + 1, A 2 + 1, A)F 4 (a, b, c - 1, c 2 ; (1 - x^l - x 2 )), (28) 



where 



2 2A 2A 

a = -n,b= -(Ai + A 2 + 1) + n - 1, a = -— , c 2 = -r^. (29) 
A A A 



Substituting (H) to eq.@ we find 

(p(q)) oc x M (l - x)"F A (-fi, p + v + 1 + n, p + ^, 1/ + ^, x 2 , (1 - x) 2 ), (30) 

where a; = sin 2 (fc^g/2A^) as before. The Appell function F 4 (a, [3, 71, 72, x, y) satisfies the fol- 
lowing system of equations 

x(l - - y 2 ^ - 2xyZ^ + [ 7l - (a + (3 + l)x}Z' x -(a + [3+ l)yZ' y - af3Z = 0, 

(31) 

y(l - y)Z" m - x 2 Z'^ - 2xyZ% + [ l2 -{a + f3 + l)y]Z' y - (a + (3 + l)xZ' x - af3Z = 0. 
In the limit q <C N we obtain after simple manipulations 

(p(q)} oc Z(u,v)\ u=2q 2 jV=0 , (32) 

where Z(u,v) is the solution of equations 

= fj,(fj,-l)Z 

(33) 

2nZl - 2^ - 2( W - v)Zl + 3Z; + k\Z = 0. 

For it 3> 1 in the first order variables u and t> are separated each from other. We have as a 
result 

(p(q))<xn + -^Z 5 (2k F q), (34) 
where Z$ is the Bessel function and 

6 = VV4 + 8M//-1). 
The asymptotic behavior of ([53]) is in agreement with general result (fZ7j). 



5 Summary and conclusion 



We have studied the asymptotic behaviour of the overlap integral between ground states of 
the BC N type Calogero Sutherland model corresponding to different strenghts of the boundary 
fields and the Friedel oscillations of the single particle density due to the presence of a boundary. 
Our results are in agreement with those obtained by applying the predictions of boundary 
conformal field theory to the finite size spectra of these systems. In particular, it has been 
established that non universal exponents showing a continuous dependence on the values of the 
boundary fields can arise in the correlation functions of boundary changing operators. 
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Figure 1: Single particle density oscillations for a system of iV = 10 particles with A = 2, // = 3 
and v — 1. 




